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Abstract 

By adapting Feynman's sum over paths method to a quantum mechanical 
system whose phase space is a torus, a new proof of the Landsberg-Schaar iden- 
tity for quadratic Gauss sums is given. In contrast to existing non-elementary 
proofs, which use infinite sums and a limiting process or contour integration, 
only finite sums are involved. The toroidal nature of the classical phase space 
leads to discrete position and momentum, and hence discrete time. The corre- 
sponding 'path integrals' are finite sums whose normalisations are derived and 
which are shown to intertwine cyclicity and discreteness to give a finite version 
of Kelvin's method of images. 



1 Introduction 

In this paper we give a new proof, using only finite sums and avoiding the need for 
analytic methods, of the Landsberg-Schaar formula for quadratic Gauss sums. The 
key idea, which springs from earlier work of one of us [Q , is to apply the quantization 
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of a quantum mechanical system whose phase space is a torus, with both position 
and momentum (and also time) discrete and cyclic. On adapting Feynman's 'sum 
over paths' method to this system a discrete and finite version of Kelvin's method of 
images is obtained. 

The formula states that for positive integers p and q 

Tv £ exp {—) = 7W) S exp \-^r) • (1) 

In number theory this formula plays a central role, underpinning key results relating 
to quadratic reciprocity and characters. The formula also promises to be a useful 
adjunct to discrete Fourier transform methods used, for instance, in some algorithms 
for quantum computing 0. In this context it essentially provides the discrete inverse 
quantum Legendre transformation from the free Hamiltonian to the free Lagrangian. 

The standard proof of the Landsberg-Schaar formula (as given, for example, in 
0) is obtained by putting r = 2iq/p + e, e > 0, in the Jacobi identity for the theta 
function: 

+oo 1 +oo 

£ e"™ V = — r e— 2/T , (2) 

n=— oo * n=— oo 

and then letting e — > 0. This method is an example of Hecke's observation that 
'exact knowledge of the behaviour of an analytic function in the neighbourhood of its 
singular points is a source of arithmetic theorems' |4| p225. 

While this method certainly establishes the truth of the formula, Hecke's insight 
notwithstanding it is somewhat unsatisfactory to have to use analytic methods and 
take limits when only finite sums are involved. The proof given in this paper does 
remain entirely in the arena of finite sums; it demonstrates the extraordinary range of 
Feynman's method of 'summing over histories', which here extends to the discrete and 
cyclic regime in a pleasing and effective way. So, in a sense, the quantum mechanics 
offers another example of Hecke's insight 

The key ingredient in this paper is quantization on a toroidal phase space. A 
first step was taken by one of us in leading to a proof of the Jacobi identity 
by using random walks to quantize a system with cylindrical phase space. The key 
idea here is the observation that the series i?(z£) = Yln=-oo e ~ 7rn lt is the trace of a 
certain quantum evolution operator exp —iHt/h, which converges in the distribution 
sense. Using path integral methods adapted to the cylinder, in particular Kelvin's 
method of images , an alternative formula for the trace is found, and equating the 
two expressions gives the Jacobi identity. This work is described in section |2| as an 
introduction to the main result in this paper. 
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With a minor change of emphasis, the formula ([[]) turns up by taking the operator 
that gives (0) (for r = it) and then 'looking at the system at discrete times r = —2i/p\ 
This idea is pursued in section |3|. In order to remove difficulties over convergence, 
and to obtain the Landsberg-Schaar formula (|I|), we change the problem to one in 
which instead of the infinite set of energy levels (0) there is only a finite number of 
energy levels E n . This is achieved, following a suggestion of Berry ||, by replacing 
the cylindrical phase space (with angular position 9 periodic but angular momentum 
L unbounded) by a torus on which both 9 and L are periodic. On quantization both 9 
and L are discrete, and moreover, as we see, discrete times are also required. Adapting 
the methods used in section ^ to prove the Jacobi identity to this new setting, again 
two expressions for the trace of the quantum evolution operator are obtained, one by 
working in the momentum basis in which the evolution operator is diagonal and the 
other by path integral methods (using a novel discrete variant on Kelvin's method of 
images); the equality of these two expressions gives the Landsberg-Schaar identity. 

This first description of the discrete and cyclic path integral method is presented 
somewhat heuristically, but in a manner which should emphasise the ideas and mo- 
tivation. In section |] full details of the toroidal quantization are given together with 
a justification of the normalisations used in section |3[ 



2 Cylindrical phase space and the Jacobi identity 

In this section earlier work of one of us |I[ is described as a prelude to the main result 
of the paper, in order to demonstrate some of the novel features of the work. 

Consider a rigid body constrained to rotate about a fixed axis (cf Schulman ||). 
Let / denote the moment of inertia and 9 denote the angle of rotation (given in 
radians so that it is taken mod 27r). An angular Schrodinger picture is used, with the 
angular momentum observable L represented as the operator 

(where h as usual stands for Planck's constant h divided by 2tt). Using the classical 
expression for the energy or Hamiltonian 

T 2 

H = TV W 

the Schrodinger equation for the wave function ip(9,t) is then 

.d^ h 2 d 2 ijj , N 

lh m = ~2iW (5) 
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Suppose that E n denotes the n th eigenvalue of the quantized Hamiltonian with energy 
eigenfunction <f> n , 

-YiW U6) = EnUel (6) 

with respect to boundary values determined by periodicity in 9 with period 2n. Evi- 
dently 

E n = MO) = e me - (7) 



The solution of fl5|) satisfying the initial condition 



+oo 

in(8-6 ) 



e 

n=— oo 



where 5 denotes the periodic delta distribution with period 2tt gives the kernel (or 
matrix element) of the evolution operator exp —iHt/h: 



1 + °° 2 

K(MA,0) = (9\e^\9 ) = ±Y. e-^e^\ (9) 



which is convergent in the distribution sense; alternatively one can give t a small 
imaginary part to restore convergence. Later we shall see that when both position 
and momentum are cyclic, with phase space a torus, the problem of convergence does 
not arise. 

From(^) we see that the trace of the evolution operator is 



1 + °° 7 2 t 

Tr (exp -iHt/h) = — ^ (10) 



with convergence in the sense described above. 

Now, as in |], we follow || adapted to the cylindrical phase space. The purpose of 
the section is to express the old idea of using the method of images and the universal 
cover of the circle in terms of path integrals, in such a way that the argument can be 
adapted to the toroidal phase space of section |3], obtaining a second expression for 
the trace of the evolution operator exp iHt/h which on comparison with (|T0|) gives 
the Jacobi identity. 

Denote by 

S(6(t)) = [ C(9,9;t)dt, (11) 
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the classical action, where C is the Lagrangian, and then we obtain the evolution 
amplitude from (9 ,0) to (9, t) as 

K{9,t;9 ,0)~ Yl e iS(e)/hbar ^ 

periodic paths from (e ,o) to (e,t) 

where '~' indicates that a normalisation factor (c.f. (|13D ) is required. In order to 
carry out the path integral sum we carry out the integration on the universal covering 
space, K, of S 1 , that is, we lift the homotopy classes to the universal covering space 
HI ||- Corresponding to paths from 9 to 9 in S 1 , we have paths from some fixed 9q 
in z/ -1 (#o) (here v denotes the covering projection from R to S 1 ) to each of the 9* 
in u~ 1 (9), where the index j runs through the fundamental group (= Z) of S 1 . It 
follows (this is of course Kelvin's method of images) that 

*■(«.**.<»= £ {iLy^im 1 - - "- 2 ^)- < 13 > 

n=— oo x ' N ' 

(The normalisation factor is taken from ||.) 

On equating the expressions for the trace in (pTOj ) and (p~3[ ) we obtain 

+00 /r> r\ - +°° 

E (ihn 2 t) I Znl \ 2 ST^ -2h 2 n 2 I 

Ti— — oo x / n=— oo 

(As already observed, care is needed over convergence, which is intended in the dis- 
tributional sense.) If we arrange that 2irl — H, r — it, then we obtain formally the 
usual Jacobi identity (0) for r d{it). 

After this work was presented |l| it was suggested by Berry || that the way 
forward to obtaining the Landsberg-Schaar formula by similar means might be to 
consider a system where the angular momentum, as well as the angle, was cyclic so 
that the phase space was compact and the quantized observables would not only be 
discrete but also have a finite range. This approach, which does indeed lead to the 
Landsberg-Schaar formula, is developed in the following section. 



3 Toroidal phase space and the Landsberg-Schaar 
formula 

In this section we adapt the methods of the preceding section to a toroidal phase 
space on which both L and 9 are periodic. On quantization this gives discreteness to 
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both these variables, and makes it possible (provided that suitable values are used 
for the various constants involved) to have only a finite number of energy eigenstates. 
Time, too, becomes quantized as emerges when evolution is considered. In this section 
we concentrate on ideas, obtaining the Landsberg-Schaar formula by proceeding in 
the manner of a physicist, while in the next section the quantization scheme is fully 
described and normalisation factors are derived. 

Suppose that the phase space of our system is a torus, with angular momentum 
L of period P and angle of rotation 9 of period 2tt. Quantization in L and 9, when 
both 9 and L are periodic, gives, respectively: 

L = nh, 6=*^, m,neZ. (15) 

The phase space is quantized to a lattice on the torus, and the small rectangles with 
sides h and ^ps fit a whole number of times into the phase space torus area. The area 
of the phase space torus is 2nP. So 2-nP = 2nNh~, (P = Nh~, N e Z) and so the small 
rectangles fit N 2 times into the phase space torus. 

There are precisely N energies, E n , 1 < n < N, and so the trace of the evolution 
operator is given by: 



1 r — iE n t 

Tr(exp — iHt/h) = — > e » 

V.7T < J 

e~ — . (16) 



n=l 
N 



1 r "V ifrn^t 



n=l 

Now we turn to the L -periodicity and show that it implies a restriction on the times, 
t, of 'looking at the system'. 

Any quantum mechanical state ip(9,t) has an expansion: 

+00 rk 2 t \ 

ip(9,t) = ^2 ake 1 ^ 6 '^) , a k = a k+N (17) 

with the usual observations about convergence. Now take the Fourier transform in 
9 of (|T7D to go from position (angle) space to momentum space. We obtain (here ip 
denotes the Fourier transform in the variable 9), 

+oo 



$(L,t) = / i;(9,t)e- vpe/h d9, 
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and we use T [e tke ] = 5(L — kh)) 

+00 

ip{L,t)= 22 a k e- — 5(L-kh), (19) 

k=— 00 

where S denotes the Dirac delta distribution, and 

+00 

<*. r — ihk t 

ip(L,t) = a k e ~5(L + P — kh) 

k=— 00 
+00 

a k e —5(L-kh). (20) 

k=— 00 

Thus periodicity of period P in L, that is ip(L + P, t) = ip(L, t), implies 

1 = e -iht{N 2 +2kN)/{2I) q\) 

so that (assuming N is even, as is required later) 

27rml , , 

t = meZ. (22) 

On substituting from in fliTf ) we obtain 

Tr(exp -iHt/h) = — J^e 



iV 



2tt 

n=0 
N 

2 it 



1 * 

-L^e-^f. (23) 



n=0 

We have to compute 

e iS ^ h , (24) 

periodic paths 

where the sum is now over all periodic paths on the torus. At the n th energy level, 
p n is constant: 

2 

E n = ^,ie n = Pn ,I9 = Pn t. (25) 



Hence, for the 'time of looking' given in fl22l 



10 <Pn- ' ' < \m\.2nl. 26 

j\h 
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On repeating the argument given in section but now restricting 9 as in (|2"E|), we 
obtain in place of (p~3|) 



m-l / T x , / . r 



= £ (55* -pU< 9 - 2 ™) 2 )> < 27 > 



27rml , , 

i = lv/r < 28 > 



n=0 

where as before 



(The normalisation factor used is taken from fll3|), and is justified in section [|; while 
it might seem simpler to derive this by the method used by Davison [[J than by the 
methods of section in fact there is a crucial step (equation (3.5) of Davison's paper) 
which would require knowledge of the Gauss sums under study.) 
It follows from (^7j) that the propagator trace is 

N \ 2 / m z 



( N \ 1 
Tr(exp -iHt/h) = {^—) $>xp 



n=0 

m— 1 



( 2h.2m-K \ 
\ Nh ) 



(^L\'Y e l ™ 2N ' m . (29) 



n=0 



So, on equating ( Pf ) and ([29]), we obtain 



JV / , T \ I m-l 



^ e-^T = f ^ " ^ e WiV/m (30) 

n=0 ^ ' n=0 

Finally in (|30D we take N = 2q, m = p, i~ l = —e m l 2 to obtain 
which is the Landsberg-Schaar formula ([[]). 



4 Normalisation of the discrete path integral 

In this section details are given of the quantum mechanical system whose path in- 
tegration is used in the previous section to prove the Landsberg-Schaar formula (|I|). 
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The necessary path integral formula is derived by adapting the original approach of 
Feynman || to the discrete and cyclic setting. Because all sums are finite in this 
case there are no convergence or other analytic difficulties, so that a precise result 
is obtained with well-defined normalisation. This confirms the validity of the more 
heuristic, but geometrically and arithmetically well-motivated, use of the Feynman 
principle in the previous section. 

The starting point is the classical phase space of the system, which is taken to be 
toroidal with coordinates 9, L of periodicity 2ir and P respectively. As Hilbert space 
for our system we take 7i = C N (where N is a positive integer) realised as 

( 2-Kin9\ 



f :R - C\f(0) = $>„exp (-^H 

n=l ^ ' 



where each a n is a complex number. 

One basis of His then plainly {f k : k = 0, . . . , N - 1} with f k (9) = ^exp(^). 
This is the angular momentum basis, that is, each f k is an eigenvector of the angular 
momentum operator L (defined now by L = —i^-§s) with eigenvalue k. Using Dirac 
notation, we write f k as \k). For simplicity we use units in which h takes the value 1. 

Another basis is {b r : r = 1, . . . , iV} with 

JV-l 



M^£exp(M^). (32 ) 



N ^ 1 V N 

k=0 v 

The inner product on 7i is defined by 

N-l 

(f\9) = J2f^*9U)- (33) 

3=0 

With this inner product both of these bases are orthonormal. 

The b r may be regarded as the position basis if one restricts the domain of the 
elements / of TC to Z N = {0, 1, . . . , iV — 1} and defines the (exponentiated angular) 
position operator x by 

9 



xf(e) = exp\2m-)f(9), (34) 
since then 

xb r = exp ^27rz— j b r (35) 

so that b r is an eigenstate of x, with eigenvalue exp (2nij^Y (In Dirac notation we 
write b r as |r).) 
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Now the proof of the Landsberg-Schaar identity essentially involves calculating 
the trace of exp —iHt (with H = |j and t = 27 ^ J as before) in these two different 
bases. From the outset we will set m — p, where p is one of the two integers in the 
Landsberg-Schaar formula. 

Method 1 is the direct way, that is, working in the momentum basis in which H 
is diagonal. 

We have Ht\k) = ^\k) so that 

Tr (exp -iHt) = £ exp (-^) (36) 
fc=0 ^ ' 



Method 2 uses discrete, cyclic 'path integrals'. We consider (r| exp — iHt\s). 
Breaking the time interval into p steps At = -, and using Sj to label basis elements 



b Si at the i step, we have 

(r| exp — iHt\s) = (r\exp—iHAt\s p -i){sp-i\exp—iHAt\s p -2) 

...(si\exp-iHAt\s) (37) 

(with summation from to N — 1 over each of the intermediate s,, i — 1, . . . ,p — 1). 
We thus need to consider exp — iHAt\si-\). Now 

(si\exp-iHAt\si-i) = (s i \k)(k\s i _ 1 ) exp ( — — J 

, q N 6XP N ) 6XP (,-lV- ) (38) 



fusing the fact that 



, . , . 1 2irkr . . . 

<r|A;> = — exp— ). (39) 



Thus 



N-l I / ni \ 

(si|exp-ii7At|sj_i) = J^-expf --(£;- (s^-s^)) 2 J 



k=0 

^2 



c W [ni [Si ). (10) 
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Again, as in section we set N = 2q, so that in particular N is even and we have 

(si\ exp -iHTlsi-i) = JT ex P i~ 2 ^ ( k ~ ( Sl ~ s ^)) 2 ) ex P f 2 " ^' 2N ~ ) 

= -Lexp^ (g --;;- l)2 ") (41) 
by equation ( p0|) of the appendix. Thus 

(r|exp-^T| S ) = (ziVr/ 2 ^... £ ^ 2^ ^,^ (42 ) 

si=0 s p _i=0 ^ ' 

where s = r, s p = s. 

At this stage we assume that p and q are coprime. Then, if p is odd, so that p 
and N = 2q are coprime, we may observe that 

/ (A;iV + /p)2\ ^ 2«t 2 / — - , . 

EE ex p 2iV ; ) = E ex p^ = (43) 

Z=0 fc=0 v F y t=0 F 



where we have used the result ( f49[) from the appendix. (Because p and iV are coprime 
the expression kN + Ip takes all Np distinct values (modA/p) as k ranges from to 
p — 1 and I ranges from to N — 1.) Also, if p = 2r is even (but p and g are still 
coprime) then 

EE- (*^) - EX>p - vTttfi)^. 

(44) 

This result allows us to insert in the expression for (r| exp —iHt\s) the extra summa- 
tions X^o 1 Sfc=o ex P {^ 7r ^ kI 2Np P ^ ) together with the compensating factor 1/y/iNp 
(if p is odd) or 1 / ^/2iNp (if p = 2r with r = 2r' even). We can overlook the case 
where p = 2r and r is odd, since in that case the Landsberg-Schaar equation is triv- 
ially satisfied. This step effectively allows the winding number k for a path to be 
shared over the p steps in the path. 

Using these facts, we see that if p and q are coprime and a = 1 when p is odd 
while a = 2 when p = 4r', 

(r| exp —iHt\s) 



11 



V f S1=0 s p _i=0 1=0 k=0 \ \i=l ^ 

V ^ Sl= o s p _i=0 Z=0 fc=0 \ \i=i \ / f 

V f Sl= o s p _i=0 Z=0 fc=0 \ \j=1 ^ 

(45) 

provided that r — s = (mod N). 

Now let Mi = si — so + I = s i — r + ^ ^2 = S2 — si + I and so on, with u p = 
s p — Sp-i + l = s — Sp-i + l. Then Y^=i u i = r — s + lp so that u p = r — s + lp — Y^ZiU- 
Thus if p and N are coprime and s\, . . . , s p _i are fixed the integer variable u p will 
take each value 0, 1, . . . , iV — 1 mod N precisely once as / ranges from to N — 1, so 
that 

(r|exp-iHT|s) 

(iJV)-*/ 2 ^ ^ Pi / /A «? \ \ / (Jfc 2 g) \ 



ui=0 « p =0 fc=0 \ \ i=l 



(M)f/ 2 ^exp ^Tri^^j by (49) 



fc=0 

Hence 



P-i 

HpN 



(46) 



AT-l 



Tr(expitft) = (s| exp iift|s) 



s=0 



= \^ PY p 2iri— — ^ 

v^iV^ V P / 

= y exp f 2**^^ (since N — 2q). (47) 
V*P£^ VP/ 
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(If p = 4r' with p and q still coprime, then the values of u p are restricted, so that 
summation over I contributes a factor exp (j^f) = V^Np. The result above is thus 
also obtained in this case.) 

Hence in all cases when p and q are coprime we have (using ([36]), 

2g-i , , 2 \ fcr- p- 1 



sr^ ( k 2 p\ J2q v-4 / (k 2 q)\ 



(4£ 



which immediately gives the Landsberg-Schaar formula when p and g are coprime. 
The general case follows on setting p = mp', q = mq' with p', q' coprime. 

Appendix 

The two formulae below may be proved by elementary means: [HI, TT] 



E( 2ni(n — s) 2 \ 
exp \ } )=y/2Fi (49) 



n=0 
2r-l 



fc=0 V 7 

Given such an elementary evaluation of the Gauss sum, together with an elementary 
proof of the quadratic reciprocity law 

S)g) -<-!>**. <«> 

where denotes the Legendre symbol and p and q are odd primes, one could 

prove (H|) by evaluating each side and then appealing to the reciprocity law, and its 
extensions (to p = 2 and p = — 1 ) . Whether such a proof would provide insight into 
why the result is true is another matter. 

Again, there is another proof of the theta function identity (g), due to Polya [|T2]1, 
which depends on an identity involving binomial coefficients, identities which are in 
turn related to the Markov chain approach to diffusion processes considered in [[[J. 

The approach adopted in this paper presupposes an elementary proof of the eval- 
uation of Gauss sums in (|4"9"D and (loTf), but not of the reciprocity law. It is perhaps 



tempting therefore to regard the present approach as a substitute, in some sense, to 
the reciprocity law, but we prefer to see it as a quantum mechanical equivalent of 
Hecke's observation, in which Feynman's 'sum over histories' in discrete time replaces 
the limiting process that derives ([]]) from (Q). 



13 



References 

[1] J.V. Armitage. Gauss sums and Markov chains. Lecture given at L.M.S. Sym- 
posium on Number Theory and Dynamical Systems, York, Great Britain, 1987. 

[2] P.W. Shor. Algorithms for quantum computation: Discrete logarithms and 
factorising. In S. Goldwasser, editor, Proceedings 35th Annual Symposium on 
Foundations of Computer Science, page 124. IEEE Computer Society Press, Los 
Alimitos CA, 1994. 

[3] H. Dym and H.P. McKean. Fourier Series and Integrals. Academic Press, 1972. 

[4] E. Hecke. Vorlesungen iiber die Theorie der Algebraischen Zahlen. Akad. Verlag, 
Leipzig, 1923. English translation: Lectures on the Theory of Algebraic Numbers 
Springer Verlag, New York 1981. 

[5] M. Berry. Private communication, 1987. 

[6] L. Schulman. A path integral for spin. Phys. Rev., 176:1558-1569, 1968. 

[7] L Schiff. Quantum Mechanics. McGraw-Hill, 3 edition, 1955,1968. 

[8] R.P. Feynman and A. Hibbs. Quantum Mechanics and Path Integrals. McGraw- 
Hill, 1965. 

[9] B. Davison. On Feynman's 'integral over all paths'. Proc. Royal Society, 
A225:252-263, 1954. 

[10] E. Landau. Elementare Zahlentheorie, first section of Vorlesungen iiber Zahlen- 
theorie. S. Hirzel, Leipzig, 1927. Reprinted: Chelsea, New York 1950. 

[11] W.C. Waterhouse. The sign of the Gaussian sum. J. Number Theory, 2:363, 
1970. 

[12] G. Polya. Elementare Beweis eine Thetaformel. Sitz. der Phys. -Math. Klasse, 
158-161, Berlin, 1927. 



14 



